Intrinsic curvature determines the crinkled edges of "crenellated disks" 
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Elastic curvature constants determine many structural and functional properties of fluid mem- 
branes. Methods to measure the mean curvature modulus have proved to be robust. In contrast, 
Gaussian curvature is an intrinsic property of a surface. Thus, measuring the relevant modulus k 
in fluid membranes remains a challenging task. Inspired from colloidal "crenellated disks" observed 
in a model system composed of hard rods, we propose a concise relation between the two curvature 
moduli and the parameters associated with the free, crinkled edges. Our approach offers a straight- 
forward way to determine k of these reconfigurable membranes, where various complex topologies 
can be nanosculpted. Further, we reveal the structure and stability of the "crenellated disks." 



I. INTRODUCTION 

A significant achievement towards a continuum de- 
scription of two dimensional (2D) fluid membranes was 
made by Helfrich, when he proposed a free energy ex- 
pansion in terms of mean (extrinsic) and Gaussian (in- 
trinsic) curvatures, the two invariants of the curvature 
tensor of a surface [l|. It follows that many physical 
and functional characteristics of membranes depend on 
the associated elastic curvature constants, such as di- 
verse vesicle shapes in fluid membranes and specifically 
red blood cells 0-111, fission and fusion processes [1, [^, 
formation of caveolae and buds in lipid bilayers or cell 
walls 0, 13 , cell division, and many more. Microscopi- 
cally, these elastic curvature constants are determined by 
lateral pressure profiles across a membrane Q , which are 
obtained through theoretical models or extensive com- 
puter simulations [lol - [T3| . As for experiments, meth- 
ods to measure mean curvature modulus have become 
reliable, such as analysis of thermal undulations along 
the membrane normal 15| or vesicle aspiration and 
tube pulling techniques (l6j]. On the other hand, the 
Gaussian curvature modulus can only be measured when 
e.g. catenoidal pores nucleate in a membrane @, or by 
the characterization of the boundary curve enclosing the 
membrane surface. These measurements utilize the well- 
known Gauss-Bonnet theorem [13, [3 ; which relates the 
boundary geometry with the total Gaussian curvature of 
the area it surrounds. However, these transformations 
of membrane conformation are in general difficult to ob- 
serve and control, posing a problem for the measurement 
of the intrinsic curvature modulus. 

In this paper we develop a theoretical model which fo- 
cuses on the boundary properties of "crenellated disks" 
observed in a model system of colloidal monolayers as- 
sembled from chiral, rod-like fd viruses 0, [H. We 
simultaneously minimize the total free energy, which is 
based on the Helfrich framework, both with respect to 
the surface and the boundary deformations of a crenel- 
lated disk. Our approach reveals a simple analytical re- 



lation between the two elastic curvature moduli and the 
parameters characterizing the crinkled edges. This rela- 
tion serves as a powerful tool to indirectly measure the 
Gaussian curvature modulus of these membranes. Then 
we determine the equilibrium structure of the crenellated 
disks. Finally, we investigate their thermodynamic sta- 
bility with respect to the perfectly flat monolayers, which 
are widely observed in the global phase diagram of fd- 
viruses mixed with a depleting agent. 

The present work is organized as follows: In the next 
section we explain the hypothesized structure of the 
crenellated disks. In Sec. IIIII we present a detailed de- 
scription of our model. In Sec. IIVI we discuss our the- 
oretical results on the structure of the crenellated disks 
and the theoretical phase diagram. Concluding remarks 
are offered in the final section of the paper. 



II. STRUCTURE OF A "CRENELLATED DISK" 

Experimentally, mixtures of monodisperse rods with 
opposite chirality are prepared from two types of viruses: 
a left-handed wild-type fd virus {fd wt) and a right- 
handed mutant fd virus, dubbed fd Y21M. These mix- 
tures were observed to form a nematic phase at a certain 
mixing ratio, and a cholesteric {Ch) phase with a re- 
duced overall strength of handedness otherwise [l^]. In 
the presence of a depleting agent, rod mixtures which 
are either achiral or have a low overall chirality assem- 
ble into the crenellated disks. These unusual structures 
are essentially flat monolayers surrounded by an array 
of local, three-dimensional (3D) cusp defects (Figs, fflja)- 
(c)) [2l|]. In the membrane interior, the rods are per- 
fectly aligned parallel to the layer normal, resembling a 
single layer of a conventional chiral smectic-A (Sm-A*) 
phase [22| . Molecular twist deformations penetrating at 
the edge are either clockwise or anticlockwise, and sur- 
round the monolayer in the form of a narrow Ch band. 
The cusps at the boundaries are the termination points of 
these alternating left- and right-handed Ch bands. Note 
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FIG. 1. Schematics of the crinkled edges surrounding 
a crenellated disk, (a) Top, (b) top-front, (c) front views 
of the orientation of rods around the cusps. Adjacent cusps 
alternate above and below the monolayer plane. In (a) the 
line tension 7i and in-plane protrusion amplitude ai of each 
free boundary are also shown. The green arrow is the axis of 
reflection. Away from the cusps the bulk of the monolayer is 
planar. For illustration purposes the concentration and the 
aspect ratio of the viruses have been reduced, (d) Schematics 
of an isolated cusp. For definitions, see text. 



that phase separation between the two types of viruses 
is observed neither in the crenellated disks nor in bulk 
mixtures. Instead of forming defects on the plane, these 
Ch regions escape into a third dimension in order to re- 
duce the orientational mismatch of rods where the edges 
meet (Fig. [Ijc)). Because molecular tilt changes direc- 
tion at each of these reversal points, adjacent cusps must 
alternate above and below the plane of the otherwise flat 
disk. Towards the defect peak, each of these Ch sec- 
tions acquires a macroscopic handedness opposite to the 
direction of molecular tilt, i.e. a Ch band becomes a left- 
handed space curve when molecular twist deformations 
are right handed, and vice versa (see the yellow bands 
in Figs.[lja) and (b)). Likewise, in a previous work rods 
with given chirality were found to form twisted ribbons 
(minimal surfaces to a double helix) of opposite hand- 
edness as a result of molecular twist deformations ^23*]. 
To maintain a continuous structure between two adjacent 
cusps, there lies an axis of reflection, as shown in Fig. [1] 
(green arrows). 
In contrast to lipid monolayers and bilayers (3l| . fd 



viruses tend to build colloidal membranes with nega- 
tive Gaussian curvature [l^ [l^ . In addition to twisted 
ribbons, which occur in the limit of high chirality and 
low depletant concentration, the strong effect of nega- 
tive Gaussian curvature on the Sm-A* membranes is evi- 
denced by two further observations: (1) Disks obtain high 
edge energies by typically reaching to mesoscopic sizes in 
diameter, however, they never fold into surfaces with- 
out boundaries, such as a sphere (a surface with positive 
Gaussian curvature), to avoid these free edges, (2) achiral 
disks stretched by optical forces again form a twisted rib- 
bon [l^ . Furthermore, a characteristic of a saddle surface 
is that it exposes a longer edge than if it were flat |24| . 
which is indeed observed both in the twisted ribbons and 
crenellated disks. Thus, we hypothesize that the forma- 
tion of cusp arrays could only be stabilized with respect 
to the flat monolayers by virtue of the negative Gaus- 
sian curvature. The saddle characteristics of the surface 
further increase by the formation of in-plane protrusions 
(Fig. HJa)) at the expense of increasing interfacial energy 
along the free edges. 



III. THEORY 

Here, we assume rotational symmetry around each 
cusp. That is, we consider a surface of revolution with 
negative Gaussian curvature for an isolated defect, where 
the peak is located at its rotation axis (Fig. HKd)). The 
variation of the free energy of a generic membrane with 
respect to the surface deformation has been calculated 
within the context of chiral lipid membranes f25l - [27| . To 
account for the formation of the in-plane protrusions as 
well, we carry out a minimization both with respect to 
the deformation of the area and the free boundaries of 
the membrane. We employ the Helfrich free energy along 
with relevant interfacial terms, and ignore the variation 
of the molecular director field to simplify the analysis. 
Then the free energy of the Sm-A* membranes is given 
as follows: 




where /if = - {2H f + knc + cr . 



Jh is Helfrich free energy density of surface bending, 
where H is the mean curvature, kQ is the Gaussian cur- 
vature, and k, k are the corresponding moduli, respec- 
tively [HQ- f is the surface tension of the membrane. 
These terms are multiplied by dS*, the membrane area el- 
ement. The interfacial terms, associated with the forma- 
tion of the free edges, are multiplied by the infinitesimal 
arc length dsi (see Fig. [Ud)). Each cusp is the termi- 
nation point of two free edges with opposite handedness, 
hence the sum in front of the closed integral. That is, two 
edge profiles are allowed to be different. First two inter- 
facial terms are; the line tension 7^, and k;,, the bending 
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rigidity of the interface. Experimentally, 7^ and ki, are 
measured from the power spectra of interfacial fluctu- 
ations on the monolayer plane. The interfacial tension 
7i is controlled by the molecular chirality of the viruses, 
as well as the depletant concentration^l^]. The natural 
length scales of the system, which determine the relax- 
ation of the 3D structure into a 2D membrane, are given 
by £,i = -y/ kb /ji . The geodesic torsion Tg^i in the last term 
of Eq. dl]) is the pseudoscalar describing the macroscopic 
chirality of the alternating left- and right-handed Ch sec- 
tions. Its prefactor is hi. This term serves as a correction 
to 7i due to the 3D structure of the cusp, and vanishes 
away from it along the edge. In Sm-A* membranes, be- 
cause molecular twist is confined to a narrow region at 
the interface, we approximate —hiTg^i to an interfacial 
term by assuming a constant molecular tilt within a nar- 
row band at the edge. This is in contrast with the case 
of chiral lipid membranes, where the tilt of the chiral 
molecules is uniform throughout the membrane, leading 
to a bulk term -hTg [2j, |23, [Hi. 

We evaluate Eq. ([T|) for a generic surface of revolu- 
tion of negative intrinsic curvature without making any 
further assumptions, except that the local height of the 
surface fir) is a single- valued function. In Euclidean 
three-space a surface of revolution of radius R is repre- 
sented by the position vector 

Y(r, If) — {r cos r sincp, /(r)} , 

where r<R, —(f>i{r)<ip<TT + (f>2{r), (2) 
and (/)i(r) , 02(?') > . 

We define a set of orthogonal coordinate vectors 
{Yr , , N} on the membrane, where N is the local 
unit normal and , Yi^ are the tangent vectors to the 
surface in radial and azimuthal directions, respectively. 
The two invariants of the surface, the first- and second- 
fundamental forms, can be defined in terms of an angle 
a(r) < and its derivatives, which satisfies tana(r) = /' 
(primes denote derivatives with respect to r). The two 
principal curvatures via these invariants are obtained as 
Ki — cosaa' and K2 — sina/r (Fig. [TJd)). Then, the 
mean and Gaussian curvatures are calculated by the rela- 
tions H = (ki + ^2)72 and kg — kiK2- The area element 
in terms of a is given by dS = rdrd4>/ cos a. 

The surface of revolution is terminated by two space 
curves confined to the membrane, corresponding to 
two free edges meeting at the peak. Each individual 
curve is represented by the position vector Yi(s) — 
Yi(r(s), 0i(s)) on the surface, where s is the coordinate 
along the arc length of the curve The curvature of 
these curves are obtained in terms of a, Pi{r), and their 
derivatives, where /3i(r) < is the angle between the tan- 
gent vector Ys_i to the i-th curve and Y^ (Fig.[T](d)). For 
a, and 0i the relation tan /3i = r cos acjj'^ always holds. 
This relation serves as an additional constraint to be 
added to Eq. ([T]) with an associated Lagrange multiplier 
\i (r) , arising from the choice of parametrization. The to- 
tal curvature squared of the free edges is given by j = 



Kg j -I- K^ i, where Kg^i = i cos a (r sin/?;)' is the geodesic 
curvature and Hn,i = kjCOS^ /3i -I- K2 sin^ jSi is the normal 
curvature of the lines [la] • The geodesic torsion of each 
curve is found as Tg i = sin/J^ cos/3i(K2 — The line 
elements in Eq. ([ij are given by dsi = dr / {cos a cos Pi) 
(Fig.lHd)). 

We perform a mean-field theory by simultaneously 
minimizing Eq. ([1]) with respect to the deformation of 
the surface and the free edges. This is done by analyzing 
the variations of the free energy ^ as a boundary-value 

problem, where gj = {a, /?2, 0i, 02} [29|]- Far from 
the rotation axis the minimum of the free energy ~ 
is achieved using the following boundary conditions: The 
surface becomes a fiat membrane where f{r) and all its 
derivatives vanish. Hence, ki, K2, and Kn,i become zero. 
Furthermore, due to two torque-free boundary conditions 
the curvature K^^i of the free edges ideally vanish, result- 
ing in a straight line in the flat membrane limit {Kg^i — 0, 
f3i = in this limit), and Xi{r) = [2S|. Plugging in 
these boundary conditions to the set of we arrive at 
the relation 

k + k = 71^1 = 72A2 , (3) 

where Ai is the maximum in-plane protrusion amplitude 
of the z-th edge (that is, '^i\ji^^ = see Fig. [TJa)). 
Eq. ([3]) is a natural consequence of the Gauss-Bonnet 
theorem, which states that the total Gaussian curva- 
ture of an area is equivalent to the total amount of the 
geodesic curvature Kg ^ and jump angles along its bound- 
aries [13, [3] ■ This ensures a practical determination of 
the Gaussian curvature modulus k, provided that mea- 
suring the line tensions 7i, the mean curvature modulus 
k, and the amplitudes Ai is relatively easy in experi- 
ments. Furthermore, it implies an elegant connection 
between Ai and A2 , specified solely by 71 and 72 , or vice 
versa. 



IV. SHAPE AND STABILITY OF 
CRENELLATED DISKS 

We apply Eq. ^ to the array of cusps surrounding 
Sm-^* monolayers composed of virus mixtures in achi- 
ral and slightly chiral regimes. In the latter, one of the 
Ch bands {e.g. associated with 72) has opposite chi- 
rality to the overall handedness of the mixture. Thus, it 
chooses an erroneous tilt direction, whereas the other Ch 
band (with 71) tilts the right way. Eventually, 72 > 71, 
and the edge of incompatible chirality becomes shorter 
(Fig. [IJa)). This results in a long-range attractive force 
between two cusps. On the other hand, when two defects 
are sufficiently close, a short-range repulsive interaction 
arises primarily from the edge curvature energy. Hence, 
frustration emerges between line and curvature energies 
of the free edges, and one should expect a finite equi- 
librium distance do{= 2Rq) between two adjacent peaks. 
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Here i?o is the equilibrium radius of a single defect. For 
convenience, we denote the edge with lower line tension 
7i, and the edge with higher line tension 72 in the fol- 
lowing. Other edge-specific parameters also follow this 
convention. 

In the membranes composed of a single type of virus, 
the dependence of the line tension on the chirality is re- 
ported in Ref. [l9| . In our current theoretical analysis, 
we model the effect of the overall mixture chirality by 
6 = 72/71, where S £ [1, 1.5]. Furthermore we set h = 0, 
since the geodesic torsion has a negligible effect both on 
the overall shape and the free energy of the crenellated 
disks. Kb is taken to be ~ lOO/csT/iw [l^, and its depen- 
dence on the mixture chirality is ignored. The lower line 
tension is measured as 71 '~ bOOksT / fim p^ . In real- 
ity, all experimental parameters depend on the chirality 
of the mixture, which will be published elsewhere along 
with complementary theoretical results [HI. 

With respect to a perfectly flat monolayer, the forma- 
tion of cusp arrays is expected to be stabilized by negative 
Gaussian curvature. Since Sm-A* monolayers have free 
boundaries and no spontaneous curvature, the Laplace 
pressure on the membrane should vanish. Thus, we con- 
sider a minimal surface {H = 0). The only nonplanar 
minimal surface of rotation is catenoid, where f(r) sat- 
isfies 



ecosh 



fir) -fie) 



(4) 



as given in Ref. [30] • Here e = —i? sin a(i?) is the 
neck radius of the catenoid (r > e), and /(e) = 

~'^lQg( j; [ i+cosa(fl)l ) is the cusp height (Fig. Elja)). We 
choose a{R) to be as small as —10^^ to suppress any non- 
analyticity at R on the surface. The boundary conditions 
at the cusp are a = — 7r/2 and /3i = While evalu- 
ating Eq. ([!]) , we ignore the effect of variations SF/S(j)i in 
the regime at ^ R. Namely, the structure of an isolated 
cusp defect is mainly governed by SF/SPi and fir). 

The phase diagram for the crenellated disks (gray- 
shaded area) and flat monolayers as functions of k and 5 
is displayed in Fig. HJa) . To obtain the phase diagram we 
calculated the difference AF between the free energy of a 
single defect (Eq. ([T])) and the free energy Ffiat of a fiat 
monolayer of an equal area S, where Ffiat = crS + 2jiR. 
Both for achiral (6 = 1) and chiral mixtures {S > 1), we 
find a first-order phase transition between the crenellated 
disks and flat membranes, as k is varied. Rq for two adja- 
cent cusps is finite at the phase boundary (see Fig.jSfc)). 
In reality, because 7; and Kb are different for fd wt-rich 
and fd Y21M-rich mixtures, the phase diagram should 
exhibit asymmetry between the left- and right-handed 
samples. However, since we model the effect of chirality 
solely by S, Fig. ^a) does not distinguish between the 
mixtures of opposite handedness, k at the phase tran- 
sition increases when Rq decreases. Essentially, if Rq 
is small, then there is less energy gain from the Gaus- 
sian curvature due to lower (see Fig. |3Ib)). Thus, 
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FIG. 2. The thermodynamic stability of the crenel- 
lated disks, (a) For a minimal surface, phase diagram for 
the crenellated disks (gray-shaded region) and flat monolay- 
ers in k versus S = 72/71. The black full curve is the first 
order phase boundary. The achiral region is at 5 = 1. (b) For 
k — 125kBT, the force curves (= \—dAF/dr\) as a function of 
r. (c) For k — 125kBT, AF as a function of r. When the min- 
ima (red dashed curve) of the curves are below AF = (black 
dotted line), then the crenellated disks are stable with respect 
to the flat monolayers. Otherwise, they are metastable. In (b) 
and (c), from bottom to top, curves are drawn for increasing 
5 with 0.1 increments, starting from achiral crenellated disks 
(3 = 1). 
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FIG. 3. The profiles of the crenellated disks, (a) The 

height /(r) as a function of r. The cusp height /(e) grows with 
increasing R. For k — 125fcsr the in-plane protrusion profiles 
are given in cartesian coordinates y versus x, corresponding to 
(b) 5 — 1.4 (chiral) and (c) 5 = 1 (achiral). The full curves are 
equihbrium in-plane protrusions, whereas the black dashed 
curve in (b) is a profile of a non-equilibrium crenellated disk. 
The asymmetry due to 5 7^ 1 of the edges is obvious in (b). Rq 
is long for achiral mixtures in (c). Eq. holds for sufiiciently 
long edges. 
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the Gaussian curvature leads to an additional repulsive 
force at low distances. The widest region of stability is 
found in achiral rod mixtures, as also observed in exper- 
iments. In the flat membrane phase shown in Fig. [2{a), 
the crenellated disks are always found to be metastable, 
in qualitative agreement with experiments. 

In the absence of the mean curvature it is straightfor- 
ward that Eq. ([3]) simplifies to 

k = 71 Ai = 72^2 . (5) 

When 5 is sufficiently low, Rq is big enough that the equi- 
librium protrusion amplitude reaches Ai, satisfying 
Eq. (O (see Fig. [3jc)). We find good qualitative agree- 
ment between the orders of magnitude of Ai and the ex- 
perimental measurements [2l|. Note that fc is a material 
parameter and is expected to be a function of the mix- 
ture composition. The measurements of k using Eq. ([5]) 
will be reported elsewhere. The value of k estimated 
in Ref. [2j] for the Sm-A* membranes was k ^ lOksT, 
which was too low to explain the case of no spherical vesi- 
cal formation. In fact k for fd virus samples must also 
be much bigger than the moduli measured for some lipid 
monolayers p3ll - l33| (again, |fc| ~ lOksT for lipid mem- 
branes). A geometry with a positive Gaussian curvature 
always possesses a positive mean curvature, which natu- 
rally holds for these lipid monolayers. Hence, in principle 
with a lower |fc| compared to the Sm-A* monolayers, lipid 
monolayers could stabilize their observed structures. 

The strength of the attractive force due to 72 depends 
primarily on 6. This is evidenced by the in-plane protru- 
sion profiles in Figs. |ni[b)-(c), as well as the force curves 
in Fig. IHb) . When S is high, there is a clear asymmetry 
between two protrusions (ai > a2), and the attractive 
force due to 72 is strong. On the other hand, when 5 
gets lower, the attractive force becomes weaker as 71 and 
72 become comparable to each other, resulting in a bigger 
i?o (see Fig. |31[c)). In the achiral regime {S = 1) a very 
weak force on the order of 10~''7i survives, since lower 
also reduces the edge length between close cusps. In this 
case, the distance between all adjacent cusps should be 
equal, as 71 =72. As stated above, due to its repulsive 
nature at low R, the negative Gaussian curvature forces 
Rq to be bigger. A simpler model ignoring the 3D struc- 



ture of the defects and accounting only for the line and 
curvature energies of the edges consistently resulted in 
smaller Rq values than the results presented here. This 
2D model only gives metastable states. Thereby we de- 
duce that the Gaussian curvature not only stabilizes the 
crenellated disks, but also ensures an accurate determi- 
nation of their structural properties. 

V. CONCLUSION 

Crenellated disks are monolayers of aligned rods which 
are decorated with an array of cusps alternating above 
and below the plane of the disk. At these cusps the 
molecular twist changes handedness. Previously, Tu and 
Pelcovits utilized a phenomenological elastic theory of 
the crenellated disks to calculate the director pattern 
of the rods on a prescribed surface with variable geo- 
metric parameters (3^ . Our work, rather, puts empha- 
sis on the calculation of the detailed optimum geometry 
which minimizes the Helfrich free energy. In our analysis 
we have ignored the variation of the molecular director 
field. The effect of the director deformations, however, 
is still accounted for by the interfacial terms in Eq. ([IJ, 
albeit crudely. The constants Kb and 7^, which are mea- 
sured in experiments, are directly related with the liquid 
crystalline distortions at the edge of the monolayer jl9| . 
Therefore, the theoretical framework presented here cap- 
tures the qualitative essence of the experimental results 
while maintaining conciseness. The primary outcome of 
this model applied to a single, isolated cusp is the connec- 
tion between the elastic curvature moduli and the quan- 
tities specifying the free boundaries, as given in Eqs. ([3]) 
and (O. In particular, the relation k = jiAi provides a 
novel recipe for the indirect determination of the Gaus- 
sian curvature modulus k of the rod mixtures in question. 
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A. Total free energy 

For a single surface of revolution terminated by two 
space curves, in terms of the angles a, Pi, and 4>i, geodesic 
curvatures Kg^i and normal curvatures K,n,i of the free 
edges, as well as the principal curvatures ki and K2 of the 
surface (defined in the main text), the total free energy 
is given by 



F = F„+F, 



H + ^edgei + Fedge2 + F\i + Px 



(1) 



where 



F, 



H 
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- dr 
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COS a cos Pi 
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-hi [K2 - Kl) 



(3) 



cos a 



= / drXi ( rcosa-^^ t&n Pi 
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The index of the edges isi = {l,2},eis the cut-off length, 
which is given by e = —Rsina{R) for the catenoid. 



B. Euler-Lagrange equations 

The variations jp, where gj — {a, Pi, P2, (pi: 4'2}, are 
equivalent to the Euler-Lagrange (EL) equations. For 
one independent and several dependent variables, the EL 
equations are given by [l| 



(5) 



dgj dx dg'j 

where / — f(r,gj,g'j) is the free energy density given 

by _F = J drf{r,gj,gj). Primes denote derivatives 
The set of EL Equations in Eq. ^ contains ten coupled 
nonlinear first-order differential equations. Evaluating 
Eq. ([U, these equations are calculated as 
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cosacos^ P1U2 ~ {k + k) h sin/3i(7i — 2ki,k^ i) 
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(12) 

In Eq. ([9]), the two constraints rcosa^^ = tan Pi are 
used. Based on Eqs. ([S]) and p^ . we conclude that A = 
Ai = A2. Then, the nine unknowns to be solved are 
a,a',^i,/3i,/32,/32,'/'i,'/'2, A. 

Note that we neglect the contribution of the variations 
1^ to the surface and edge deformations, as explained 
in the main text. Therefore, the structure in our approx- 
imation is governed solely by a, a', /3i, /3^, /32, /32- Iii the 
case of the minimal surface, since a and a' are calcu- 
lated analytically (see main text), the equations govern- 
ing /3i,/3[ and /32,/?2 decoupled from each other. 



C. Boundary conditions 

The torque-free boundary conditions are given by 
Eqs. 0, and dH) being equal to zero 0. Suffi- 
ciently away from the rotation axis of the surface, in 2D 
membrane limit, a — 0, a' — 0, and a" — 0. Then, 
Kn.i = 0. Plugging in these results to Eq. ([9|), we find 
that sin Pi = Ai/r. This rather expected result applies 
to a straight line in polar coordinates, where the shortest 
distance between the origin and the straight line is Ai 
(or, the in-plane protrusion amplitude). When hi = Q as 
in the achiral limit, then the same relation still holds, and 
this time the torque-free condition = is satisfied. 
The geodesic curvature Kg^i of a straight line vanishes in 
2D, hence ^7 = 0. In the 2D membrane limit A must be 

zero. This is best seen when a curve with the free energy 
given in Eq. ([3]) is minimized on a plane. The resulting 
EL equations are independent of (pi, then A = 0. Hence 
= 0, away from the rotation axis of the nonplanar 

surface. Plugging in these results into Eqs. pH)) and pT|) . 
we find that k+k = 71^1 = 72^2. For a minimal surface, 
k does not contribute to this equation. 

At the cusp, ideally, the boundary conditions are a — 
— 7r/2 and Pi — Q (again, we ignore dependence of 

the structure). In the case of the minimal surface, since 
e « lO^^i?, we take a — 0.95 which nearly corresponds 
to a 70° slope. This slope is steep enough as evidenced 
in Fig. 3(a) in the main text. 
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